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ABSTRACT 

A power-law time-dependent lightcurve for active galactic nuclei (AGNs) is expected by the self-regulated 
black hole growth scenario, in which the feedback of AGNs expels gas and shut down accretion. This is 
also supported by the observed power-law Eddington ratio distribution of AGNs. At high redshifts, the AGN 
life timescale is comparable with (or even shorter than) the age of the universe, which set a constraint on the 
minimal Eddington ratio for AGNs on the assumption of a power-law AGN lightcurve. The black hole mass 
function (BHMF) of AGN relics is calculated by integrating the continuity equation of massive black hole 
number density on the assumption of the growth of massive black holes being dominated by mass accretion 
with a power-law Eddington ratio distribution for AGNs. The derived BHMF of AGN relics at z = can fit the 
measured local mass function of the massive black holes in galaxies quite well, provided the radiative efficiency 
^0.1 and a suitable power-law index for the Eddington ratio distribution are adopted. In our calculations of the 
black hole evolution, the duty cycle of AGN should be less than unity, which requires the quasar life timescale 
tq > 0.5 giga-years. 

Subject headings: accretion, accretion disks — ^black hole physics — galaxies: active — quasars: general 



1. INTRODUCTION 

There is evidence that most nearby galaxies contain massive 
black holes at their centers, and the central massive black hole 
mass is found to be tightly corr elated with the velocity disper- 
sion of the galaxy dFerrarese & Merritt 2000; iGebhardt et alj 
I2OOO1) . or the luminosity of the spheroid component of its hos t 
galaxy (e.g., iMagorrian etanil998t iMarconi & Huntll2003h . 
These correlations of the black hole mass with velocity dis- 
persion/host galaxy luminosity were widely used to estimate 
black hole masses, and to derive the mass functions of the 
central massive black holes in galaxies (e.g. .lYu & Tremaine 
2002; Marconi et al. 2004; Shankai- et al. 20041; iTamura et al. 
2006. .Graham et al.. ,2007i: .Shankar, ,2009) . On the other 



hand, quasars are powered by accretion onto massive black 
holes, and the growth of massive black holes could be domi- 
nantly governed by mass accretion in quasars. The massive 
black holes are the refore the active galactic nuclei (AGN) 
relics ( ISoltan|[l982l) . and the luminosity functions (LF) of 
AGNs provide important clues on the growth of massive 
black holes. The black hole mass function (BHMF) of AGN 
relics can be calculated by integrating the continuity equa- 
tion of massive black hole number density on the assump- 
tion of the growth of massive black holes being dominated 
by mass accretion, in which the acti vity of massive black 
holes is described by the AGN LF (e.g.,'Cavaliere et al.l'1971l; 
Soltan 1982; Chokshi & Turner 1992; Small & Blandford 
19921: IShenI 120091: iLi et al.1 120091) . There ai-e two free pa- 
rameters: the radiative efficiency ?7rad and the mean Ed- 
dington ratio A = iboi/i'Edd for AGNs, adopted in most of 
such calculations on the cosm ological evol ution of massive 
black holes (e g..lYu & Trema ine 2002; Mar coni et al.ll2004t 
IShankaretalJl2004l) . The derived BHMF of AGN relics in 
this way is required to match the measured local BHMF at 
redshift z = by tuning the values of two parameters rj^^i 
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and A, which usually requires almost all AGNs to be accret- 
ing close to the Eddington limit (e.g..lYu & Tremaind 120021 ; 
IMarconi et al]l2004t[Shankar et alJl2004lK 

In principle, the mean Eddington ratio for AGNs A is not 
a free parameter, which can be estimated from a sample of 
AGNs with measured black hole masses. One of the most 
effective approaches for measuring masses of bl ack holes in 
AGNs is the reve rberation mapping method (Petersonlll993l ; 
iKaspi et a^l2000^ . Using the tight correlation between the 
size of the broad-line region and the optical luminosity es- 
tablished with the reverberation mapping method for a sam- 
ple of AGNs, the black hole masses of AGNs can be easily 
estimated from their optical luminosity and width of broad 
emission line. The mean Eddington ratio A ~ 0.1 at z 
0.2 to ~ 0.4 at z 2 was derived from a large sample of 
AGNs with the analyses of the Slo an Digital Sky Survey 
(SDSS) by McLure & DunloD (2004) (also see 'Warne r eTal] 
2004; Kelly et al. 2010). Kollmeier et al. (2006) pointed that 
the samples selected from SDSS are heavily weighted to- 
ward hig h-luminosity objects d ue to the limited sensitivity 
of SDSS. lKoUmeier et alJ (l2()06h estimated the Eddington ra- 
tios of AGNs discovered in the AGN and Galaxy Evolution 
Survey (AGES), which is more sensitive than the SDSS. The 
derived Eddington ratio distribution at fixed luminosity is 
well described b y a single log-nor mal distribution pe aked at 
- 0.25 (also see iShen et a l. 2008; Trump et alj|2009l) . How- 
ever, some other investigations showed that the Eddington 
ratios of loc al AGNs spread over several orders of mag- 
nitude (e.g., Ifol I2OOI iHopkins et all l2006h . The Edding- 
ton ratio distribution for AGNs exhibits a power-law distri- 
buti on with an exponential cutoff at a high Eddington ra- 
tio dMerloni & Heinzl 120081; IHopkins & Hernquistll2009l) . or 
a power -law distribution with a n additional log-normal com- 
ponent (iKauffmann & Heckmani i2009). Such a power-law 
Eddington ratio distribution is qualitatively consistent with 
the self-regulated black hole growth scenario, in which the 
feedback of AGN expe l s gas and shut down accretio n (e.g., 
IHopkins et all l2005aibt IHopkins & HernquistI l2009h . This 



2 



means that an AGN with bolometric luminosity Lboi may con- 
tain a relatively small black hole accreting at a high rate or 
a more massive black hole accreting at a lower rate. In this 
work, we adopt a power-law Eddington ratio distribution with 
an exponential cutoff at a high ratio to derive the BHMF 
of AGN with a LF, with which the continuity equation for 
black hole number density is integrated to calculate the cos- 
mological evolution of BHMF of AGN relics. The resultant 
BHMF of AGN relics is constrained by the measured local 
BHMF. The conventional cosmological parameters Om = 0.3, 
r^A = 0.7, and Hq = 70 km s~' Mpc"' have been adopted in 
this work. 

2. THE EDDINGTON RATIO DISTRIBUTION FOR ACTIVE 
GALACTIC NUCLEI 



iHopkins & Hernquisj (I2009h suggested that the quasar 
lightcurve can be described by 



dt 



dlogX 



A 



^peak 



exp 



A 



^peak 



(1) 



where A = Lboi/^Edd (^Edd = 1.3 x lO^'^Mbh/M© ergs s '), tq 
is the quasar life timescale, and the parameter Apeak describes 
the peak luminosity of quasars. This power-law lightcurve is 
consistent with the self-regulated black hole growth model, 
in which feedback produces a self-regulating ^fdecay^s or 
^blowoutqs phase after the AGN reaches some peak lumi- 
nosity and begins to expel gas and shut down accretion (e.g., 
IHopkins et al.ll20()5al) . This lightcurve can be translated to an 
observed Eddington ratio distribution C(A), 



C(A): 
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where C/ is the normalization, if the switch-on of AGN ac- 
tivity is balanced with switch-off of AGN activity and tq 
is significantly shorter than the age of the universe at red- 
shift z- The power-law Eddington ratio distribution is consis- 
tent with those derived with samples of near by AGNs (e.g., 
Heckman et all l2004t lYu et alj gOOS; Hop kins & HernquistI 



2009h . Such an Eddington ratio distribution has a lower cutoff 
at A = Ajnin.o, below which the sources are no longer regarded 
as AGNs. In this work, we adopt A = Amin.o = 10""* in all the 
calculations. 

At high redshifts, the quasar life timescale tq is comparable 
with (or even shorter than) the age of the universe at redshift 
Z. In this case, the time after the birth of the first quasars is 
so short that most of them are still very luminous (i.e., accret- 
ing at high rates), and the lower limit on the Eddington ra- 
tios for AGNs should be significantly higher than Amin.o. The 
first quasars are predicted to have formed a t Zfq ^ 10 (e.g., 
iHaiman & LoebllioOUlBromm & Loebl2003l) . with which we 
can estimate the minimal Eddington ratio X'^^^ for AGNs at 
redshift z as 
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Apeak 
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d\ogX = t{z), (3) 



where f(z) is the age of the universe at z measured from 
Zfq = 10 when the first quasars formed. For simplicity, we 
adopt Amin = max[An,in,o, A^;,,(z)] in all our calculations on the 
black hole evolution. In Fig. [T] we plot the minimal Edding- 
ton ratios as functions of redshift z for different quasar life 
timescales tq. 




Fig. 1 . — The minimal Eddington ratios as functions of redshift z for dif- 
ferent AGN life timescales: tq = 10^ (red), 2.5 X 10* (green), 5 X 10* (blue), 
and 7.5 X 10* years (black), where A^jip o = lO"* is adopted. 

For a given black hole mass function A^AGN(A^bh,z) and Ed- 
dington ratio distribution, the AGN LF $(z,Lboi) can be cal- 
culated with 



NagMu , z) ^1^1^ ((X)d log A Mpc-3(logLboi)-' , 
c/logLboi 



(4) 

where M^u/Mq = Lboi/(l-3 x 10^*^ A), A^AGN(Mbh,z) is the 
AGN BHMF at z, and the Eddington ratio distribution C(A) 
is given by Eq. (O. In this work, we assume that the AGN 
black hole mass function has the same form as the AGN LF, 



A^AGN(Mbh,z) = 



A^AGN,o(z) 



(Mbh/M*b)A»' +(Mbh/M,\)A»2 



Mpc ^(logMbh) 
(5) 

where the parameters A^agn.o, -^bh' Ani and /3m2 are to be de- 
termined. Substituting Eq. (|5]l into Eq. (|4|l, we can calculate 
the AGN LF with a given Eddington ratio distribution C(A) 
provided the values of the four parameters in the AGN BHMF 
are specified. In this work, we tune the values of these param- 
eters till the observed LF can be well reproduc ed by that cal- 
culated with Eq. (|4|. We adopt the LF given bv lHopkins et alj 
( l2007bl) . which is calculated by using a large set of observed 
quasar luminosity functions in various wavebands, fr om the 
IR thro ugh optical, soft and hard X-rays (see Hopkins et alj 
I2007bl for the details). 



$(Lboi,z) = 



(Lbol/i*)'^'+(ibol/i,)'^^ 



Mpc-^(logLboi)"', (6) 



with normalization 0*, break luminosity L*, faint-end slope 
7i, and bright-end slope 72. The break luminosity L^, evolves 
with redshift is given by 

logL, =(logL,)o + ^L,lC + ^L.2?^ + ^L,3^^ (7) 

and the two slopes 71 and 72 evolves with redshift as 



and 
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71 = (7i)o(-; ) 

1+Zref 
2(72)0 



(-_l+Z_\K7T.l 1 ( 1+Z ',<t^^.2 



(8) 



(9) 
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The parameter ^ is 



C = log 



1+z 

1+Zref 



(10) 



and Zref = 2 is fixed. The other parameters of this LF are 
as follows: log 0, (Mpc"^) = -4.825 ± 0.060, [logL,(3.9 x 
lO^^erg s-i)]o = 13.036 ± 0.043, k, , = 0.632 ± 0.077, k, 
-11.76±0.38, 



'<-L3 



-14.25 ±0.80, 



(7i)o 



: 0.417 ±0.055, 



^71 



-0.623 ± 0. 132, (72)0 = 2.174 ± 0.055^J3^j = L460± 
0.096, and k^,,2 = -0.793 ±0.057 (see Hop kins et alJl2007U 
for the details). This LF includes both the Compton-thin and 
thick sources. 

We give some examples of the AGN BHMFs derived from 
an AGN LF with a given Eddington ratio distribution at dif- 
ferent redshits z in Fig. |2l It is found that the AGN LF can 
be well reproduced by the calculations from the AGN BHMF 
with the form given in Eq. (|5]l. 
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Fig. 2. — Top: The comparison between the best-fitted LFs (lines) cal- 
culated form the AGN BHMF with a given Edd ington ratio distribution 
and the LFs (dots) given by HopkinFetalJ )200 7b) at different redshifts: 
z = (red), 3 (green), and 5 (blue). The parameters: f}/ = 0.3, Xpcnk = 2.5, 
and TQ = 5 X 10* years, are adopted. Bottom: The derived AGN BHMFs at 
different redshits: z = (red), 3 (green), and 5 (blue). 



3. THE COSMOLOGICAL EVOLUTION OF MASSIVE BLACK HOLES 

The evol ution of massive black ho le number density is de- 
scribed by (iSmall & Blandfordl[T992h 



dn(Mkh,t) d 



dt 



dMu 



■[n(Mbh,?)<M(Mbh,?)] = 0, (11) 



where the black hole mass function «(Mbh,f) is in units of 
Mpc~^Mj^iJ, < M(Mbh, > is the mean mass accretion rate for 
the black holes with Mbh, and the efforts of mergers of black 
holes are neg lected (see Shankar et al. 2009; Cao & Li 2008] 
for the detailed discussion). The black hole mass function 
A^(Mbh,0 used in this work is in units of Mpc~"'(logMbh)~\ 
Equation ( fTTT ) can be re-written as 



1 dN(M^^,t) d 



Mbh dt 



N{M^,k,t)<M(M^,k,t) > 



Mbh 



= 0, 
(12) 

since N{M\,h,t) = In 10Mbh«(Mbh, 0- 

It was suggested that the standard thin accretion disk will 
transit to an advection dominated accretion flow (ADAF) 
when the dimensionless mass accretion rate m is lower 
than a critical value merit (tn = M/Medd; Mem = 1.3 x 
lO^'^Mbh/O.lMgc^) dNarayan& Yilll995h . The radiative effi- 
ciency for ADAFs is significantly lower than that for standard 
thin disks, and it decreases with decreasing mass accretion 
rate m. It was suggested that the radiative efficiency rj^-id can 
be described with 



'7iad,0 



if ffi > nir 



cnt5 



(13) 



where mem = 0.01 is adopted in this work (see lNaravanll2002l 
for the detailed discussion and the r eferences there i n), an d 
the parameter s = 1 is suggested in iNaravan & Yil (Il995h . 
The calculations of the ADAFs surrounding rotating black 
holes in the general relativistic frame showed that the value 
of s is in the range of ^ 0.2- 1.1 depend ing on the value 
of black hole spin parameter a (Xu & Cao 2010). In all our 
calculations, we adopt s = 0.5 (e.g., Merloni & Heinz 2008]; 
[Draper & Ballantvnell201C)b . The mean mass accretion rate 
for the black holes with Mbh at redshift z can be calculated 
with 



N(Mbh,z)<M(Mbh,z)> 

f C(A)A^AGN(Mbh,z)(Mbh/M0)ALEdd,0(l - r?rad) 
J Vr^dC^ 



J log A, 



(14) 

where NACN(Mbh,z) is the AGN BHMF at z, and the radiative 
efficiency rji-^ is given by Eq. ( fT3] l. 
The black hole evolution equation can be rewritten as 



1 9A^(Mbh,z) dt d 



Mbh dz 



dz dMbh 



N{Mbh.t) < M{Mbh,t) > 



Mbh 



(15) 

As described in §2, the AGN BHMF A^AGN(A^bh,z) can be cal- 
culated with a given Eddington ratio distribution (|2]i and the 
AGN LF (|6]l, and the mean mass accretion rate can be calcu- 
lated with Eq. ( fT4b . Integrating Eq. ( fT5] l over z from z^ax with 
derived AGN BHMF A^AGN(Mbh,z) and suitable initial con- 
ditions at Zmax, the cosmological evolution of massive black 
holes is available. The duty cycle 6 of AGN is defined as 



S(Mbh,z) = 



NACN(Mbh,Z) 

N(Mbh,z) ' 



(16) 



which is required to be less than unity. In all our calcula- 
tions, we assume that the duty cycle i5 = 0.5 at Zmax = 5. There 
are three free parameters, ryrad.o, A, and Apeak, in our calcu- 
lations, which are tuned to let the BHMF of the AGN relics 
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at z = fit the measured local BHMF given in lShankar et alj 
( l2009h . This local BHMF encompasses the range of several 
estimates of BHMF with different methods ("Bell et al."2003| 
Marconi et al. 2004; Shankar et al. 2004; Ferrarese & Foq 
2005t iTundo et al.ll2007i: iHopkins et al.i ,2007a). We find that 
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the final results are insensitive to the initial conditions at Zma, 
because the fraction of local black hole mass accreted at very 
high redshifts can be neglected. 

The results for the evolution of the total BHMFs and AGN 
BHMFs with redshift z are plotted in Fig. |4]for different val- 
ues of the model parameters. In Fig. |5j the duty cycles 5 of 
AGNs are plotted as functions of black hole mass Mbh at dif- 
ferent redshifts z- In this work, the AGNs accreting at > merit 
are referred as bright AGNs, in which radiative efficiently 
accretion disks are present. We plot the radiative efficiency 
evolving with redshift in Fig. |6] 
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Fig. 3. — The local BHMFs aX z = 0. The shaded area encompasses several 
estimates of the local BHMF (see Fig. 5 in Shankar et al. 2009. and refer- 
ences therein). The lines represent the BHMF of the AGN rehcs calculated 
with AGN life timescale tq = 5 X 10** years. The red solid Hnes are the re- 



sults with , 



; 0.3 and A, 



peak 



: 2.5, while the red dashed and dash-dotted 



lines represent the results with Xp^ak = 1 and 3.5, respectively (/?/ = 0.3). The 
blue dashed and dash-dotted Hnes represent the results with /3; = 0.2 and 0.4, 
respectively (Ap^^k = 2.5). 



4. DISCUSSION 

Our estimates of the lower limits on the Eddington ratios 
of AGNs show that Amin increases with redshift z, which im- 
plies that the mean Eddington ratio is higher at high redshifts. 
At high redshifts, the quasar life timescale tq is comparable 
with (or even shorter than) the age of the universe at red- 
shift z, and most of the AGNs are therefore still very lumi- 
nous (i.e., accreting at high rates). The mean Eddington ra- 
tios for these AGNs are relatively higher than those at low 
redshifts. The estimates of the Eddington ratios for AGNs 
show that the mean Eddington ratio increa ses with z (e.g., 
iMcLure & Dunlopll2004 iWarner et aI]|2004D . which is quah- 
tatively consistent with our results. 

Unhke most of the previous works (e.g.. | Yu & Tremaind 
I200I iMarconi et aUlIOOl iShankar et alll2004D . in which a 
single mean Eddington ratio is adopted as a free parameter, 
we use an Eddington ratio distribution in our calculations. 
Such a power-law Eddington ratio distribution for AGNs 



is expected by the s elf-regulated black hole growth model 
dSilk & Ree s 1998; Hopkins et al.ll2005iJ) . which is also sup- 
ported by the Eddington ratio estimates for AGNs (e.g_^ 
Heckman et al. 2004; Yu et al. 2005; Hopkins & Hernq uisl 
20091: Ikauffmann & Heckmanl2009l) . 
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Fig. 4.— The BHMFs at different redshifts: z = (red), 1 (green), 2 (blue), 
3 (yellow), and 4 (black). The solid lines represent the total BHMFs, while 
the dashed lines are for AGN BHMFs. The parameters for the Eddington 
ratio distribution: /3; = 0.3 and Apeai^ = 2.5, are adopted. The shaded area is 
the measured local BHMF. In the upper panel, the results are deiived with 
AGN life timescale tq = 2.5 X 10** years, while the results with tq = 5 X 10** 
and 7.5 X 10** years are plotted in the middle and lower panels, respectively. 

There are three free parameters, T^i-ad.o. and Apeak^ in our 
calculations on the evolution of massive black holes. The 
resultant local BHMF for the AGN relics is sensitive to the 
value of the peak luminosity of AGNs. It is found that the 
measured local BHMF can be well reproduced by the BHMF 
of the AGN relics at z = calculated in this work, if the 
three parameters: ?7rad = O.Il, /?/ = 0.3, and Aneak = 2.5, 
are adopted (see Fig. [Jl. IHopkins & HernquistI mm sug- 
gested that Pi ~ 0.3-0.8 based on their self-regulated black 
hole growth model calculations, and our calculations also pro- 
vide a useful constraint on the value of /?;. Our results show 
that the peak Eddington ratio of AGNs ^2.5 is required for 
modeling the local BHMF, which implies that a small frac- 
tion of AGNs are accreting at slightly super-Eddington rates. 
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This is consist ent with the estim a tes for different samples 
of AGNs (e g.. iWarner et al.l l2004t IWull2009t lAi et alJl20Tot 
iWiUott et aljl20T0lK 
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Fig. 5. — The AGN duty cycles as functions of black hole mass Mbh at 
different redshifts: z = O(red), l(green), 2(blue), 3(yellow), 4(black), and 
5(cyan). The soHd lines represent the duty cycles of all AGNs (i.e., the 
sources with A > A„,in), while the dashed lines represent the duty cycles of 
bright AGNs (i.e., the sources accreting with m > merit)- The parameters 



for the Eddington ratio distribution: 
calculations. 



Pi = 0.3, Apeak = 2.5, are adopted in the 



The AGN BHMFs can be calculated from the AGN LF ei- 
ther with a distribution of Eddington ratios or a single mean 
Eddington ratio. The derived BHMF of AGN relics at redshift 
z = is required to match the measured local BHMF, which 
usually leads to almost all AGNs be accreting close to the Ed- 
dington limit in the calculations with a single mean Eddington 
ratio. This is in consistent with the observations of AGN sam- 
ples (e;gjH^002; McLure & Dunlop 2004; Warner et al. 
2004t iHopkinset al.,.2006: .Kollmeieret al..,2006: .Kellv et ah 
20101) . The main difference between this work and the pre- 



vious works is that a power-law Eddington ratio distribution 
instead of a single mean Eddington ratio is adopted in this 
work. The AGN BHMFs derived in this work are larger than 
those calculated with a single Eddingt on ratio close to unity 
adopted in those previous works (e . g.. lYu & Tremaind 120021 : 
iMarconi etal] I2004t iShankar et 5] 120041) especiallv at low 
redshifts, because the mean ratios for the power-law Edding- 



ton ratio distributions used in this work are significantly lower 
than unity. We find that the AGN BHMFs are larger than the 
total BHMFs for relative small black holes at low redshifts 
if the quasar life timescale tq < 5 x 10^ years (see Fig. |4|, 
i.e., the duty cycle for AGNs S > 1 (see Fig. |5]), which is un- 
physical. These results imply that the quasar life timescale 
Tq > 5 X 10^ years, which is qualitatively consistent with 
the results derived either from theoretical calcula tions or ob- 
serva t ional d ata (Di Matteo et al. 2005; Hopkins & Hernquisi 
l2009HKellv"et al . 2010). As an Eddington ratio distribution is 
adopted in our calculations, the AGNs with black hole mass 
Mbh may be accreting at very low rates. We also calculate 
the cosmological evolution of the duty cycle 6h for the bright 
AGNs accreting at m > merit as functions of Mbh (see Fig. |5]). 
We find that the bright AGN duty cycles d\, are significantly 
higher than 5\, only at low redshifts, while they converge at 
high redshifts (z > 2), which is caused by the value of Amhi in- 
creasing with redshift z (see Fig. [TJ. It is found that the bright 
AGN duty cycles are higher at high redshifts, and they de- 
crease with increasing black hole mass at low redshifts, which 
are qualitatively consistent with those obtained in the previous 
works (e.g., Mai'coni et al.l l2004l; iMerlonil 12004; Wang et alj 
I2OO6; Shankar et al. 2009). Our results show that the massive 
black holes were grown earlier than their less massive coun- 
terparts (see Fig. |4|, which is similar to that found by Merloni 
(|2004). The black holes with Mbh < lO^M© were dominantly 
grown at redshifts z < 1 (see Fig. |4|. 

The Eddington ratios derived with brigh t AGN samples 
usual l y exhibit a log- n ormal distribution dKoUmeier et alj 
l2006t iShen etaH llOOa IKellv et al.1 l20Tol) . while those de- 
rived with the sample s containing fainter sources show a 
power-law d istribution (Heck man et al.ll2004 lYu et alj|2005l; 
iHopkins &' Hernquist 2009), or a power-law distribution with 
an additional lo g-normal com ponent at high Eddington ra- 
tios (K auffmann & Heckman|[2009 ). The peak and the dis- 
persion of the log-normal Eddington ratio distribution were 
found to be al most independen t of the black hole mass and 
redshift (Kollm eier et alj|2006l) . or they have very weak de- 
pendence on the black hole mass and redshift for a large 
AGN sample sele ct ed from the Sloa n Digital Sky Survey 
dShen et alj 120081) . ICao & Lil ( l2008h calculated the evolu- 
tion of massive black holes us ing a log-normal Ed dington ra- 
tio distribution derived by Kol lmeier et alj (120061) . and found 
that the measured local BHMF always cannot be reproduced 
by their model calculations unless a black hole mass depen- 
dent radiative efficiency is assumed (see iCao & Lil 12008 , for 
the details). The sources in a large sample of nearby galax- 
ies are separated into two populations by the age of galax- 
ies, which show a power-law and log-normal Eddington ra- 
tio distributions respectively (Kauffmann & Heckman 2009). 
The power-law Eddington ratio distribution for the subsam- 
ple exhibits a similar slo pe, -0.8, which is indepen dent 
of black hole mass (see iKauffmann & HeckmanI 120091 for 
the details). In this work, a power-law Eddington ratio dis- 
tribution independent of black hole mass and redshift is as- 
sumed in our calculations, which seems to be a reasonable 
assumption. We note that the slope of the power-law Ed- 
dington ratio distribution derived in this work is flatter than 
that of the distribution derived with a subs ample of nearby 
galaxies in IKauffmann & HeckmanI (120091) . The observed 
Eddington ratio distribution of all sources in their sample is 
a mixture of a power-law and log-normal distributions, the 
slope of which in the range of A - lO""*- 10-2 becomes flat- 
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ter than that of the po wer-law distribution d erived with the 
subsample (see Fig. 5 in lKauffmann & Heckm an 2009) (note 
logL[OIII]/Mbh ^1.7 corresponding to the Eddington rate). 
In this work, we use a power-law Eddington ratio distribution 
with an exponential cutoff at a high ratio in all our calcula- 
tions in order to avoid inducing additional parameters, which 
can simul ate the observed Eddingto n ratio distributions quite 
well (see ^ Hopkins & Hernquisll I20b9i, for the details). This 
distribution can also describe the main feature of the observed 
Eddington ratio distribution for the whole sample given in 
iKauffmann & Heckman (2009). The investigation on the evo- 
lution of massive black holes by adopting more realistic Ed- 
dington ratio distribution (e.g., a power-law+log-normal Ed- 
dington ratio distribution) will be carried out in our future 
work. 
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tion rate m is considered. We perform calculations with differ- 
ent values ofs (from s = 0, i.e., a m-independent 77iad, to 5 = 1), 
and find that the final results are insensitive to the value of s. 
It is found that the radiative efficiency evolves little with red- 
shift (see Fig. |6]). This is because the fraction of mass accreted 
in ADAF phases is small compared with that in b right AGN 
phase (Cao 2005; Hopkins et al. 2006; Cao 2007t IXu&Caol 
120101) . iWatarai et al.. (,200 0)ai''s calculations on the slim discs 
showed that the radiative efficiency will not deviate signifi- 
cantly from that for standard thin discs if Lboi/i-Edd 2, which 
implies that the present adop ted radiative efficiency indepen- 
dent of Eddington ratio (fT3b is indeed a good assumption for 
the sources with m > m„h- 

The Eddington ratio distribution is derived from a power- 
law quasar lightcurve in this work. The situation becomes 
complicated at high redshifts when the switch-on of AGN ac- 
tivity is not balanced with switch-off of AGN activity, and the 
Eddington ratio distribution (|2| may not be valid for the AGNs 
at high redshifts. However, the black holes were dominantly 
grown up at z < 3 (see Fig. |4]i when a power-law Eddington 
ratio distribution is believed to be a good approximation. It 
has therefore not affected much on our final results. 

The comparison of the calculated BHMF of the AGN relics 
at z = with the local BHMF gives a lower limit on the quasar 
life timescale tq. The BHMFs of the AGN relics at z = are 
insensitive to the adopted value of tq, while they become sig- 
nificantly different at relatively high redshifts (e.g., z = 1, see 
Fig. m. They can be compared with the measured BHMF at z, 
which will set a further constraint on the quasar life timescale 
Tq, provided the BHMFs are available at relatively high red- 
shift z. This is beyond the scope of this work. 



■ 1 2 3 4 5 

z 

Fig. 6. — The radiative efficiencies as functions of redstiift z- The lines are 
derived with different quasar life timescales, tq = 2.5 X 10** (green), 5 X 10* 
(blue), and 7.5 X 10** years (black), respectively. The solid hues are the results 
calculated with .v = 0.5, while the dashed lines are for s = 1. 

In this work, we adopt a m-dependent radiative efficiency 
77rad, in which the accretion mode transition with mass accre- 
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